Enhancement of thermal photon production in event-by-event hydrodynamics 
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Thermal photon emission is widely believed to reflect properties of the earliest, hottest evolution 
stage of the medium created in ultra-relativistic heavy-ion collisions. Previous computations of 
photon emission have been carried out using a hydrodynamical medium description with smooth, 
averaged initial conditions. Recently, more sophisticated hydrodynamical models which calculate 
observables by averaging over many evolutions with event-by-event fluctuating initial conditions (IC) 
have been developed. Given their direct connection to the early time dynamics, thermal photon 
emission appears an ideal observable to probe fluctuations in the medium initial state. In this 
work, we demonstrate that including fluctuations in the IC may lead to an enhancement of the 
thermal photon yield of about a factor of 2 in the region 2 < pT < 4 GeV/c (where thermal photon 
production dominates the direct photon yield) compared to a scenario using smooth, averaged IC. 
Consequently, a much better agreement with PHENIX data is found. This can be understood 
in terms of the strong temperature dependence of thermal photon production, translating into a 
sensitivity to the presence of 'hotspots' in an event and thus establishing thermal photons as a 
suitable probe to characterize IC fluctuations. 



I. INTRODUCTION 

The recent observations of large elliptic flow [l|, Hand 
energy loss of high energy partons in the medium 0] 
strongly suggest the formation of Quark-Gluon Plasma 
(QGP), a hot and dense state of quarks and gluons, in 
ultra-relativistic heavy ion collisions. The results from 
the Relativistic Heavy Ion Collider (RHIC) during the 
last decade as well as the recent results from the Large 
Hadron Collider (LHC) @ are consistent with theoreti- 
cal calculations which assume early thermalization and a 
nearly ideal fluid behaviour of the system. 

The evolution of the system from initial formation to 
final breakup is reflected in various different observables. 
Electromagnetic radiation has long been proposed as the 
most promising and efficient probe to characterize the ini- 
tial state of the hot and dense matter produced in the col- 
lisions 0. Unlike bulk hadrons, which are emitted from 
the freeze-out hyper-surface, photons come out from ev- 
ery phase of the expanding fireball, suffer negligible final 
state interactions (since electromagnetic coupling a <^ 
strong coupling as), and carry undistorted information 
about the medium conditions at their production points. 
The thermal emission of photons from the evolving sys- 
tem shows a very strong temperature dependence: the 
high pt photons are emitted mostly from the hot and 
dense early stages of matter when the hydrodynamic flow 
is weak, whereas those at lower px are emitted from the 
relatively cold later parts of the system where a large 
build-up of radial flow boosts their transverse momen- 
tum. Thermal photons with px > I GeV/c, are there- 
fore expected to provide a glimpse of the early part of 
the expansion history when the fireball is in the plasma 
phase, complementary to the spectra of bulk hadrons 
which reflect the late conditions close to decoupling. 

Initially direct photon spectra were studied mostly to 
extract the temperature of the system. In recent times, 
several other photon observables have shown significantly 



more potential to understand the dynamics of the system, 
for example, intensity interferometry of photons for space 
time information of the evolving system [T], |8j , momen- 
tum anisotropy of the initial partons as well as formation 
time of QGP from elliptic flow of thermal photons [3, |K| , 
study of jet-quenching thr ough photons produced from 
jet medium interaction |lll4l3l | etc. In addition to that, 
a higher level of sophistication has been reached in the 
estimation of prompt photons from next to leading or- 
der (NLO) perturbative QCD calculations [ij where the 
result explains the direct photon spectrum frorn p+p col- 
lisions at RHIC quite wefl for pr > 1 GeV/c Also 
photons from jet-plasma interaction, investigated in de- 
tail recently have been suggested as an important 
contribution to the direct photon yield at > 2 GeV/c. 

Ideal hydrodynamics with a smooth initial density pro- 
flle has been used successfully to model the evolution 
of the system at RHIC where the experimental data for 
hadron spectra and elliptic flow are well reproduced up 
to Pt ^ 2 GeV/c (see e.g. [H, However, for any 

given event there are inhomogeneities in the initial en- 
ergy density profile and therefore event-by-event fluctu- 
ating density profiles are more realistic than a smooth 
initial density distribution. Thus, the important ques- 
tion is whether these event-by-event fluctuating initial 
conditions (IC) have any sizable effect on the physical ob- 
servables if they are computed by averaging over many 
final states rather than considering an averaged initial 
density profile. 

It has been shown recently that hydrodynamics with 
event-by-event fluctuating IC reproduces the experimen- 
tal charged particle elliptic flow even for the most central 
collisions at RHIC [l^ which was underestimated by all 
earlier hydrodynamic calculations using smooth IC. Sev- 
eral other studies have shown that fluctuating IC give a 
better agreement with the experimental charged particle 
spectra at hig h px by increasing the number of particles 
there [3 Il9j | These fluctuations may also help to un- 
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derstand the various structures observed in two-particle 
correlations [13, [2l| . Since thermal photons are sensitive 
to the initial temperature, they are especially suitable for 
probing fluctuations in the IC. 

Thermal radiation is predicted to be the dominant 
source of direct photons in the range 1 < pT < 3 
GeV/c [H, [12] ■ Photon results using smooth profile 
and latest rates [1^ fall, however, well below the ex- 
perimental data points (see Fig. 4 of Ref. [l^) in this 
Pt range. In the present work we report a substantial 
Pt dependent enhancement in the production of ther- 
mal photons for pT > 1 GeV/c with an event-by-event 
fluctuating IC relative to a smooth initial-state averaged 
profile in the ideal hydrodynamic calculation. Conse- 
quently, our results with fluctuating density profile show 
a better agreement with the PHENIX experimental data 
for 2 < < 4 GeV/c. However it is to be noted that the 
main aim of this paper is to investigate the effect of initial 
state fluctuations on the production of thermal photons 
rather than explaining the experimental direct photon 
spectrum, which contains significant contributions from 
other sources apart from thermal radiation. Also im- 
portantly, we find that the enhancement scales with the 
initial fluctuation size. 



II. EVENT-BY-EVENT HYDRODYNAMICS 
AND INITIAL DENSITY PROFILE 

We use the event-by-event hydrodynamical model de- 
veloped in [l^, which has been applied successfully to 
study the elliptic flow with fluctuating initial states. The 
standard two-parameter Woods-Saxon nuclear density 
profile is used to randomly distribute the nucleons in the 
colliding nuclei. Two nucleons i and j from different nu- 
clei are then assumed to collide whenever they satisfy the 
relation 



(a;, - -I- {y, - y^f < IHIL^ 



(1) 



where crjvAf = 42 mb is the inelastic nucleon-nucleon cross 
section at ^Jsnn = 200 GeV, and Xi , yi denote the posi- 
tions of the nucleons in the transverse plane. 

The initial density profile is taken to be proportional 
to the number of wounded nucleons (WN), where entropy 
density s (or the energy density e) is distributed in the 
{x, y) plane around the wounded nucleons using a 2D 
Gaussian as a size function, 
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We refer to this as sWN (eWN) profile. The parameter 
if is a fixed overall normalization constant and <7 is a free 
parameter determining the size of the fluctuations. 

The effective interaction radius for the colliding nu- 
cleons, ^tTjvAr/7r/2 ~ 0.6 fm, sets a natural order of 
magnitude for the size parameter. Similarly to Ref. [l8| . 



the default value of a is taken as 0.4 fm, however we test 
the sensitivity of the results to a by varying a from 0.4 
to 1.0 fm. We consider the 0-20% most central Au-l-Au 
collisions where the value of Nwn fluctuates from 391 to 
197 according to the Monte Carlo Glauber model. This 
corresponds to an average impact parameter of about 4.4 
fm. 

Using this procedure we thus find intra-event fluctua- 
tions (i.e. 'hotspots' and 'holes') in the IC in a given 
event, as well as inter-event fluctuations in the total 
entropy between any two events. We call events with 
larger than average entropy production "hot events" and 
those with smaller than average entropy production " cold 
events" . 

For Au-|-Au collisions at ^/sWn = 200 GeV the con- 
stant K is taken as 102.1 fm^^ (37.8 GeV/fm) for the 
sWN (eWN) proflle, which reproduces the initial total 
entropy of Ref. [t?] when averaging over many initial 
states in central collisions for cr = 0.4 fm. As in [Tsj . 
the initial time is flxed at tq = 0.17 fm/c motivated by 
the EKRT minijet saturation model [l3, [11] for all the 
events. 

Assuming longitudinal boost invariance, we solve the 
(2-1-1) dimensional ideal hydrodynamic equations using 
the SHASTA algorithm ^ which is able to prop- 
agate shock waves possible with the irregular IC. The 
Equation of State (EOS) which shows a sharp crossover 
transition from plasma phase to hadronic matter is taken 
from [13 . The temperature for freeze-out is taken as 160 
MeV, which reproduces the measured px spectrum of pi- 
ous [2^ well for both the profiles eWN and sWN. 



III. THERMAL PHOTONS 

The rate of thermal photon production in the QGP 
depends on the momentum distribution of the partons 
which are governed by the thermodynamical conditions 
of the matter. 

The quark-gluon Compton scattering and quark-anti- 
quark annihilation are the leading order processes for 
photon production in the partonic phase, whereas in next 
to leading order, bremsstrahlung process contributes sig- 
nificantly to the production [2^ . In a hot hadronic gas at 
a temperature of the order of pion mass, tt and p mesons 
contribute dominantly to the photon production because 
of the low mass of pious and the large spin iso-spin degen- 
eracy of p mesons, making them the most easily accessible 
particles in the medium [301 • The leading photon produc- 
ing channels involving n and p mesons are tttt — >■ p7, 
TTp — >■ 7r7, and p — )■ tttt"/. 

We use the rates R = EdN/Spd'^x of [H for the 
plasma and those of [HI for the hadronic matter which 
at present can be considered as the state of the art. The 
transition from the plasma rates to the hadronic rates is 
assumed to happen at a temperature of 170 MeV in this 
study. The total thermal emission from the quark and 
hadronic matter phases are obtained by integrating the 
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FIG. 1: (Color online) Thermal photons from a smooth (long 
dashed blue line) and fluctuating (solid red line) sWN initial 
density profiles along with experimental direct photon data 
from PHENIX [H, ^ for the 0-20% centrality bin. The 
fluctuation size parameter a — 0.4 fm here. Also shown are 
the hadronic contributions for the smooth IC cases. 



emission rates over the space-time history of the fireball 
as 

EdN/(fp = j d'^xR{E*{x),T{x)y (3) 

where E*{x) — p'^Up_{x). The 4-momentum of the pho- 
ton is — {pt coshY, pt cos (I),pt sin (j),pT smhY), and 
the 4- velocity of the flow field is u'^ = jT{coshr],Vx,Vy, 
sinh?;) with = v'?j.=v1+Vy. The vol- 

ume element is df^x dr dx dy drj, where r — (i^— z^)^/^ 
is the longitudinal proper time and r]= tanh~^(z/t) is 
the space-time rapidity. The photon momentum is 
parametrized by its rapidity Y, transverse momentum 
Pt, and azimuthal emission angle (j). 



IV. RESULTS 

We compare our results from smooth and fluctuating 
sWN initial density profiles with PHENIX data for 0- 
20% centrality bin in Fig. [TJ The smooth initial density 
profile here is obtained by taking an average of 1000 fluc- 
tuating initial profiles which is enough to remove essen- 
tially all fluctuations. To test whether there is a residual 
dependence on the IC averaging procedure we vary the 
value of a from 0.4 to 1.0 fm. The thermal photon spec- 
trum from the resulting smooth initial density profiles is 
almost independent of a in the low px region. Only at 
high pt (> 4 GeV/c), a change in the value of a from 0.4 
to 1.0 decreases the photon production there by 20-30%. 

Photons from the fluctuating IC scenario are estimated 
by averaging photon spectra from sufficiently large num- 
ber of random events. We make sure that addition of 
another very hot or very cold event does not change the 
results significantly. 
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FIG. 2: (Color online) Ratio of photon production from fluc- 
tuating and smooth IC at different pr as a function of the 
fluctuation size parameter a. 



We find that in the entire px range the spectra are 
dominated by the QGP radiation both for the smooth 
as well as for the fluctuating IC. This is not due to our 
choice of large freeze-out temperature T/ (160 MeV): We 
have checked that a prolonged hydro evolution towards a 
smaller value like Tf — 120 MeV contributes significantly 
below 1 . 5 GeV / c but is not substantial for pT > 2 GeV/c, 
as shown by the brown dashed dotted line in Fig[T] Pho- 
tons from (only) the hadronic phase are shown by dotted 
and short dashed gray lines for Tf values of 160 and 120 
MeV respectively for the smooth IC. Both these lines fall 
well below the plasma and the total contribution. It is to 
be noted that a bag model EOS produces more hadronic 
photons near the transition region compared to a lattice 
based EOS, however the contribution from the hadronic 
phase is not significant for pT > 2 GeV/c even for a bag 
model EOS. We also note that the Bag model yield of 
photons at high px is about 50 % smaller than the re- 
sults shown in here. 

The slope of the photon spectrum from the fluctuating 
IC is about 10% flatter compared to the slope from the 
smooth one in the region 2 < px < 4: GeV/c and the 
two results differ almost by a factor of 2 (see Fig. [2]) for 
Pt > 2 GeV/c. In the region < 2 GeV/c, the fluctuat- 
ing case produces 20-40% more photons than the smooth 
proflle. The thermal emission of photons is exponential 
in temperature and linear in radiating volume. As a re- 
sult, the hotspots in the fluctuating events produce more 
high Pt photons compared to the smooth profile. On the 
other hand, the low px part of the spectrum comes from 
the relatively cold more volume-dominated later plasma 
stage, where the emission is not affected significantly due 
to the change in the IC. Thus the difference between the 
two profiles is small near ~ 1 GeV/c and enhanced 
towards higher px- 

It may be noted that a similar enhancement of the 
thermal photon production by a factor ~ 2 in the re- 
gion 1 < Px < 4: GeV/c using the same smooth IC can 
also be obtained by choosing a smaller initial time tq. 
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In this case, one needs to reduce the the value of tq to 
half of its original value • A systematic study of ther- 
mal photons at different centrality bins has the potential 
to distinguish between the effects of fluctuations in the 
initial state and the formation time. 

Similar to the photon results, hardening of the spectra 
is observed also for hadrons where the fluctuation-driven 
radial flow makes their spectra harder compared to the 
spectra obtained from a smooth IC. This is because the 
hadrons are emitted from the much later stage of the sys- 
tem expansion, i.e. the surface of freeze-out [18, 32J. For 
photons, the difference between the fluctuating and the 
smooth IC is an early time effect (due to the presence 
of hotspots) when the radial flow is still very small. Al- 
though the photons from the late hadronic phase are af- 
fected significantly due to the increased radial flow in the 
fluctuating IC, the sum spectrum, which is dominated by 
plasma radiation in the entire pt range, remains almost 
unaffected. 

Thus, the production of thermal photons shows a very 
strong dependence on the averaging procedure. We 
observe that the results for fluctuating IC show bet- 
ter agreement with the experimental data compared to 
earlier studies with smooth initial profile while leaving 
enough space for the prompt contribution and jet con- 
version photons. This may also affect the photon elliptic 
fiow results [oj . 

The sensitivity of the results to the fluctuation size 
parameter a is shown in Fig. [2l where the ratio of the 
results from fluctuating and smooth IC are plotted as a 
function of a for different values of px ■ The enhancement 
is maximal for higher values of px and lower values of 
a which clearly shows the dominance of the hotspots in 
high Pt photon production. For cr > 0.8 fm, profiles with 
fluctuating IC do not show pronounced differences to the 
smooth case as the ratio approaches 1. 

We also notice that the photon results are quite sensi- 
tive to the initial profiles; the results for sWN are harder 
than those with the eWN. The entropy initialization leads 
to more pronounced maxima in the energy density dis- 
tribution (causing hotter initial state and harder spec- 
tra) than their energy initialization counterparts which 
was observed in some earlier studies also 33]. However, 
the photon spectrum from fluctuating IC considering the 
eWN profile shows a similar enhancement relative to the 
smooth profile as observed for sWN in Fig. [TJ 

In Fig. [3] the time evolution of the photon results for 
smooth (solid lines) and fluctuating (dashed lines) IC are 
shown at px values of 1, 3, and 5 GeV/c, to illustrate 
the dynamics leading to the results presented in Fig. [T] 
For the fluctuating IC we choose two very hot events 
(red and orange dashed lines) and two relatively cold 
events (blue and green dotted lines) whereas the smooth 
IC event is shown by the solid black lines in the figure. 
The dN/dpTdrdY results clearly demonstrate that most 
of the high px photons are emitted quite early where the 
photon emission rate drops by two orders of magnitude 
within a time period of about 2 fm/c for pt > 3 GeV/c. 
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FIG. 3: (Color onIine)Time evolution of thermal photons for 
transverse momentum values of (a) 1 GeV/c, (b) 3 GeV/c, 
and (c) 5 GeV/c. Results are compared with an initial state 
averaged event and four different random events. 



The cold events refiect a relatively smaller system life- 
time. The most interesting observation is that the cold 
events produce more photons at > 3 GeV/c com- 
pared to the smooth profile due to the presence of a few 
hotspots. 



V. SUMMARY AND CONCLUSIONS 

In conclusion, we find a large enhancement in the pro- 
duction of thermal photons for event-by-event fluctuating 
IC compared to a smooth profile in the ideal hydrody- 
namic calculation. For px < 2 GeV/c, the production is 
increased by 20-40% for the fluctuating initial conditions, 
whereas for px > 2 GeV/c, it increases by about a factor 
of 2, which gives better agreement with PHENIX data 
in the region 2 < px < 4: GeV/c. A small value of the 
plasma formation time (0.17 fm/c) ensures that a fraction 
of the pre-equilibrium photons is accounted for which is 
an important source of direct photons at large values of 
transverse momentum and hard to estimate separately. 
The contributions from the initial hard processes and jet- 
matter interactions are significant for px > 2 GeV/c, and 
should be added with the thermal radiation in order to 
explain the data in the entire px range shown in Fig. [T] 
This enhancement is a result of the strong temperature 
dependence of the thermal photon emission rates which 
specifically probe the hotspots. We checked that the pho- 
ton spectra for the sWN profile are harder than the spec- 
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tra for eWN profile, however the quantitative difference 
between the fluctuating and smooth IC remains similar 
for both the cases. A more detailed investigation of the 
wounded nucleons and binary collision profiles for both 
entropy and energy initializations combined with photon 
elliptic flow from the fluctuating IC would be very valu- 
able. The difference between the smooth and fluctuating 
IC is expected to increase for peripheral collisions as well 
as for lower collision energies. We postpone these issues 
for a future study. 
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